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7^^ irest=istmicro+iblksmicro<!:: L^ wlu (irest : lenmacro, istmicro: istmicro 
+ iblksmicro-l) ^ l-^TCg-C^ U y Km-eJ^^53'tiJ L/:: ^ O ^'n y i> tB 

-oL-T^ K®<7)fIi^wlumicrotl=? — i"^o lenmicro=(lenm 

aro-irest+numthrd)/numthrdh wlumicro(lenmicro+iblksmicro, iblks 
micro) n tf— L^ lenblksmicro=lenmicro+iblksmicrO(!: "f'^o ^LX^ T^J- y 
•7°S6 3-t:\ -If 7';V-^>LUmicro^lI^^/aii-o ^titliSV^Tfi. wlumicrod inm 
icro+iblksmicro, iblksmicro) ^^ttiS't"© -^T" y 7°S 6 4 trti^ wlumicrotC^ 
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[0 0 6 0] 

;^-r-7 7°S6 l<7)*lJg|f56^~NO<7)t^'a-t;(i. ;^xy7°S 6 5tCi3V>-r. nwidthmi 
cro>=3* i b 1 ksm i cromax itz\t^ nw i dthm i cro<=2* i b 1 ksm i cromax^S^S ^ lllf i" -5) 
o X-r -yys 6 5 0^i3J8ff;{i^YE S<7)^i'^tC{i. ^ x -y 7" S 6 6 tCisv^T. nwidth 
micro2=nwidthmicro/2> istmicro2=istmicro+nwidthinicro2> nwidthmicro3= 
nw i dthm icro-nwi dthm icro2i: L. 7.7- y -fS 6 8 t^iltfo ^^t" y 6 5 OflJIff 
:^>^'N0<7)^'^tC(i, ^'.x y ys 6 7 i3i5V>T. nwidthmicro2=nwidthmicro/3, is 
tmicro2=istmicro+nwidthmicro2, nwidthmicro3=nwidthmicro-nwidthmicro2i L 
. y yS 6 8 tiattfo X-r ^ T'S 6 8 istimicro(i. -eo^^, nwidt 
hmicro t L Tnwidthmicro2 ^i!t L T -9- y'JV — ^ y interLU 'k^-Xf)^-f o 
[0 0 6 1 ) 

7.7- -j-fS^ 9tCi5V^T{i. wlu(istmicro:istmacro+nwidthmicro-l)<7)gi5^?: 
Mff1-^o itL{±. — o<7);^ix y K-eMif-rtL{f^:^-e*^o itLlCwlu(istmic 
ro: istmacro+nwidthmicro2-l, istmicro: istmacro+nwidthmicro2-l) OTH:^tT^!j 
<Oittf?iJ«r;£;6-ib^D/;^><7)-eMifi-'2>o x 7 0 tCioV>T(i. wlu(istm 

icro2: lenmacro, istmicro2: istmicro2+nwidthmicro3-l) ^wlu(istmicro2: 1 
enmacro, istmicro: istmicro2-l) Xwlu(istmmicro: istmacro+nwidthmicro2-l, ist 
macro+ 

nwidthmicro2: istmacro+nwidthmicro-1) ^5|V^'CM§f't"'5>o ^(Dh^^ l-^TCB 
trXlx y YWLXi^m^^n\^Xmmn-t^<. ^.y-y-f^ 7 1 t2;Mv^-C{±. istm 
icrot LT> istmicro2^ nwidthmicro t LTnwidthmicro3^iS L.T-^^-:/;v — ^ > 
interLU ^DfCfaiL. > ^|l-T-r'2>o 

[0 0 6 21 

mi \ 9(i. •9-r;^-^>LUmicro<7)7 0-Ti^'l)o 

y -f ^ 1 ^ \,z^\j^X ^ A (k. n/numnord) . wlul (lenblks, iblksmacro 
)^ wlumicrodeniblksmicro, iblksmicro) LT^tt-So ddt?^ wlumicr 

>XwlumicroiZ^m^tifzU6^<7)U}^m^nno istfi. LU^«|-r'2> 60 
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n y ^iliM-e^)>i)o iblksmaxti. rn -y^ft^^il (8®JS) T'*) l9 . i:fxJ^±/h^ 
< L'^v^o wlumicroti;^ l^ u/ K^t;?li§:<h LTig^tL^o 
[0 0 6 33 

^•r y T'S 7 6 tCi3V^T(±. nwidth<=iblksmaxy6-S75^*^J»f-r'2>o v^-ry^S 
7 6O^JSIf^5N0c0^{±. ;^X'y7'S 8 8 t^Btfo X-r -y-T'S 7 6<7)flJ»f56^'YE 
SOJt-^Uti, 7.7-yy°S 7 7 Kii\^^X. i=isttLT. X -f y 7°S 7 8 

i<ist+nwidth75^Si^^^*y©fi-'2>o y y°S 7 8 (Dmm^^NO(Dm^i^it^ 

[0 0 6 43 

y.-r 8 1 tc^»/^T(i. ^(7)t:4^ MigT^^^^y - Vm^-DMMy'^ y ^ 

•r -;/7°S 8 2 0*lJif75^^YE S<7)t#'^ti(i. 7" y^'S 8 3 tiisv^T. ^y-Ktr 

n#;i/cfiffi^se^ijipt;^^ ;^t- 8 6 iciitj-o y>x -y ys 8 2 tcistt 

u y ^U^hX^m^ho AtL#x./diS*K^iJipt^^&^L. ^^.x-y-T'S 8 6 
[0 0 6 53 

7.^ y-f^ 8 5-e{±. y- YWiXm.'m'LXWX^-^^ y Y^^-thJ-Yis-h"^ 
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8 6 IcioV^Tti. ^Mif L> x-r •;'7'S 8 Ttc^v^T. i ^ij i:^f TMif^:^?: 
HifL. i = i+l^LT. ;^-r 8 tcM-So 

[0 0 6 61 

>^7^ .y ■7°S 8 8 ICiJV^Tfi. nwidth>=3*iblksmax«b^ v^ii, nwidth<=2*iblksm 

7°S 8 9 tCioV^T^ nwidth=nwidth/2. ist2=ist+nwidth2 <b L> Xt^ ^T'S 9 1 
liZmtfo X-r-y-y'S 8 8 (Ommi^^N 0<Dm'^l>^ii. 7.y^yy°S 9 OtliSV^T. nw 
idth2=nwidth/3. ist2=ist+nwidth2. nwidth3=nwidth-nwidth2^ L. X-r y7°S 

9 HZMiso :^v-y'fS9 l -eti. ist(i-e<7)f nwidth^ LTnwidth2^?ll5: 
tLXr&LX, ■9-7>'W--f->LUmicro^Df LFm-To X 7" y 7° S 9 2 -Cti. wlumicr 
o(istmicro: istmacro+nwidth2-l, istmicro+nwidth2: istmicro+nwidthmicro-l) <D 
pR:$3" ?r M^'f" -2) o wlumicro(istmicro: istmacro+nwidth2-l, istmicro: istmacro+n 
width2-l)<7)TH:^^f^iJ<7)^^f^iJ^:£;&^'b*L/v:^<7)-eMif-r^o >^-r y-:^S 9 

3 "Cfi^ wlumicro (istmicro2: lenmacro, istmicro2: istmicro2+nwidthmicro3-l) 
^wlumicro(istinicro2: lenmacro, istmicro: istmicro2-l) X wlumicro (istmicro: i 
stmacro+nwidth2-l, ist+nwidth2: ist+nwidthmicro-1) $:?|V^-CMifi"-i>o X-r y 
9 4 lC;fev^T(iist^ LTist2. nwidtht LTnwidth3^^l^i: LX^ifmL 
-itaMV-^^LUmicro^O^tfai LT. -^-^^Iz-f^ > =^&lt'i>o 
[0 0 6 71 

112 Oti. •^•7')V~^yhtoc<Dyu-X$>^o 

-7 T'S 1 0 0 tCiSV^T. A (k. n/numnord) , wlul (lenblks, iblksma 
cro). buf s( lenblks, iblksunit). bufd(lenblks, iblksunit) li^-^^ttT . # 
/ — K<7) i #S <7)itiiiblksunit<7)yn y ?:numnord'f@:^;fa7t: ^) <D<DMMy^ y ^ 

^"JM^iblksmacroX iblksmacro=t «9 TOgR^^rnumnordlStC^tiJ L/i ^ <7) 

[0 0 6 81 

y ys 1 0 1 T'f±. nbase=(i-l)*iblksmacro ( i fiO^TO ^ > 
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)i'—'f(DMy)MLMWO N ibs=nbase+K ibe=nbase+iblksmacro^ len=(n-ibe)/n 
UIIlnord^ nbase2d=(i-l)*iblksunit^ ibs2d=nbase2d+K ibe2d=ibs2d+iblksunit 

^ L^ i^im'T'— -^^it^ lensend=(len+iblksmacro)*iblksunit i: "^^o 
[0 0 6 93 

;^T- y T'S 1 0 2 tlisv^r^ iy=lhLx 7 "/S 1 0 3 icioU^T. iy>nuiiino 

mi^^o ^7-yy°S 1 0 3 (Dmm^^l^ 0(Dm^l^ii. Xf-yys 1 0 4tcioV^-C 
> mmi-^M^. ^m-t^^^'^-^i^^o ■t^^t>h. idst=mod(nonord-l+iy-l,n 
umnord) +l^ isrs=mod(nonord-l+numnord-iy+l, numnord) +lt"t"'?>o -^xyT'S 1 

^*>> icp2ds=(idst-l)*iblksunit+l^ icp2de=icp2ds+iblksunit-l^ bufd(l:len 
+ iblksunit, 1: iblksunit)-^wlul (1: len+iblksmacro, icp2ds: icp2de) h-fho 1 
>k7t ^^TsUy K^-e^#J L T#y^ KT'3£^!J n tf - ^ o 
[0 0 7 03 

y^-r vT's 1 0 6-e(±. ^y- K-eiH^Mi-^o bufd<7)i*i#^idst#ioy- 

Kic^^f). bufsic^jf-r^o ^^.x-yT'S 1 0 7-Cill^«O^T^^*>> Xy-y?" 
S 1 0 8 H*5V>T. ^^'j TlWIfl^^lJ^^o y 1 0 9T-{i. KT^g^ 
tlfij t)#v^TV>^ili5iblksunit<7)^t:^-7""n -7 , ■^(DT<DU^(^y'^ y ^ . 

<D 1 ^TC a ^numnord-e:^tl| L L T^c ^ ^ (DU^ i^^^<7) y - K 

1S:S=B(7)^><^) =&7c:^*o/zfP:^H^&^-r^o A (ibs: ibe, ibs2d: ibd2d)^bufs 
(1: iblksmacro, 1: iblksunit), icps=ibe+(isrs-l)*len+K icpe=isps+len-U A( 
icps: icpe, ibs2d: ibe2d)'^bufs(iblksmacro+l : len+iblksmacro, 1 : iblksunit) t 

-r^o ki-fi 1 >x7t^^^ 1^ vm^^^ Lx^^ V y Yxnm^^wr 

^ o 

[0 0 7 13 

X-r -7 ys 1 1 0 H:joV^T(i. iy=iy+l^ LT. 7.r- y 7°S 1 0 3 KM^o 
1212 1 {±. •9-7';U-^>exrw07n--e^^o 
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[0 0 7 2] 

Xx^yS 1 1 5t3ioV>T(i> A (k. n/numnord) . wlul (lenblks. iblk 

smacro) LT^tt^o wlul (1: iblksmacro, 1: iblksmacro) Kii. LU 

^m^tif-zUnU^'k^y - V:/}mmtXW'oX^^^o nbdiag=(i-l)*iblksmacr 

'eab'2>o \^^y h<7)A:h.#x.<^'lf^7&^ ip(nbdiag+l:nbdiag+iblksmacro) 

[0 0 7 3] 

X-r -y/S 1 1 6t*{i. nbase=i*iblksunit ( i iiWmLi> t(D-iry')^-^ 
y p L\J(Oji 4 y — irows=nbase+K irowe=n/nuinnord> 

len=(irowe-irows+l)/numthrd> is=nbase+(nothrd-l)*len+K ie=min(irowe, is+ 
len-1) t-t^o >^7- yzfSl 1 7 ix=ist-t^o 
[0 0 7 4] 

X-r y7"S 1 1 ST'ti. is<=ieX-h^-^^^:^^'kmmir:ho X-r ^ T'S 1 1 8 <7) 
^|^iJBIf^^N0(7)i#'a-t^{ix X-f yy'S 1 2 Sl^zMiso X-ry^S 1 1 8 <7)^{^ijif75^^Y 
ESO^-^tKi, X7"y-7"'S 1 1 9tli5V^-r> nbdiag=(i-l)*iblksmacro. j=nbd 
ag+li:LT. :^r- yfS 1 2 0 tli3V^r> j<=nbdiag+iblksmacro-e^^^^S^^^ 
*IJifi-'2)o >^x-;'7'S 1 2 0<7)*IJ»f7!>^^NOO:fei^tI{i. >?.x-y7°S 1 2 4 tlsttf 
o Xx -y ys 1 2 0(7)it£iJif;^JYE SOJt^tKi, x ^ S 1 2 UliS v^T. ip 

S 1 2 3 tcittfo XT^yyS 1 2 1 co^flJUfT&^YE S<7)J^'^t;{i, T.^ y'fS 1 2 
2Ci3V^T. A (j, ix) tA (ip(j), ix)^AtL#x.T. -y ^ S 1 2 3 Hji 
tso Xx y 7°S 1 2 3 tlioV^Tfi. j=j+li: LT. v^t" y ys 1 2 0 t3M'i)o 
[0 0 7 5] 

;^x y/S 12 4 tlisv^rti. ix=ix+li: L. y T'S 1 1 8 
>^-r y 1 2 5 HiSV^Tfi. ^^'';T|WI^ (^/-K. 1/ y K) *^/?>o 
;^T- y7°S 1 2 6lli3V^T(i. A (nbdiag+1 inbdiag+iblksmacro, is: ie)— TRL ( 
wlul (i : iblksmacro, 1: iblksmacro))~lx A (nbdiag+1 Inbdiag+iblksmacro, is:i 

e) ^^y-K. ^^7.V"j K-esEff-T'&o c:v:t% trl (b) (i. ^t^iJB<^)TH:^ 
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oT, fettle 
[0 0 7 6] 

1112 2:S.CF[II2 3 ti^ ■9--7';U-f^>iiimcbtO-7n--e^;i>o 

X-r -y/S 1 3 0 HiSV^T, A (k^ n/numnord) ^ wlul (lenblks. iblksm 
aero) . wlu2 (lenblks, iblksmacro) LT^lt^o wluHC, -fuy^t^^ 

iblksmacrooyn y ^ ^ L U55-|!? L7t:|g:f:-e. ^t^r~n y i -g-OTfi^'n y 

[0 0 7 7] 

>^f- yys 1 3 1T'{±. nbase=(i-l)* iblksmacro ( i (iP^tFUi L ^) i: O^J--?';^ 
-^>p LUOp< ^ V^U-T'O^tJ^LIll^) . ibs=nbase+K ibe=nbase+iblks 
macro> len=(n-ibe)/numnord^ nbase2d=(i-l) *iblksunit> ibs2d=nbase2d+l^ ib 
e2d=ibs2d+iblksunit> n2d=n/numnord^ lensend=len+ iblksmacro <i: jiKMT"'" 
^ ^ (i , nw 1 en= 1 ensend* i b 1 ksmacr o t "t* •?> o 
[0 0 7 8] 

;^7^ y ^"S 1 3 2 joV^T. iy=l (^^JS^'fi^lx^) > idst=mod(nonord, numnor 
d)+l (iS'^Tfc/— (^y— F) ) ^ isrs=mod(nonord-l+numnord-l, numnor 
d)+l (^#x/-F#-^) . ibp=idsti:1->i>o 
[0 0 7 9] 

vJ.-r -y7°S 1 3 3ti43V^-r. iy>numnord-C^^ ;J>^S75^^^i^lJKff-f ;^T- -y7°S 
1 3 3<7)il^lJSlf>6^^YE ScT^^-^tCfi. ^'D^y-^y^m'iho 7.7- y -f^ \ ^ Z <D 
m^if-^0(n>^'k\iZ.{t^ ^x-^rs 1 3 4Ci3V^-C. iy=l^^S75^=^^J8ff-t-'S>o ^ 
x-yys 1 3 4 <7)flJ»f7{»^YE SO*&^t3{±, X-r ■;/ 7'S 1 3 6 ilsitfo v -f 
S 1 3 4 c7)i|£iJg|f755^NO<7)tt'g-H(±, ^x -yyS 1 3 5 l^^v^T^i^jf <7)t#^^ 
Oo ^x-y7°S 1 3 6-e(±. iy=numnord -e^^^^^^T&^^flJKff-r'^ 
o >?.7" y/S 1 3 6<7)flJ8ff;&5'YE SO^i^lifi. y 1 3 8 tcittfo 
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■y/S 1 3 6^flJif;^^'N0O^i^tC(i, ^^x y^S 1 3 7 lliSV^T. m^m^n 
wlul<7)|*l#* (^:lirn y i> 4>-^i6-r) |%<7)y-K (y-K#-^idst) Km 

[0 0 8 0] 

Xf- yT'S 1 3 8tCj3V>T. wlul<7)7"-i5^ ^j^ofcMffO-i^vv 3 >^tf^-t 
-So ibp=mod(ibp-l+numnord-l,nuinnord)+l> ncptr=nbe+(ibp-l)*len+l ( 1 -^tc 

pmmSrPft/ttit* o ^^wlul^^l §vjtt"o >^X 'y7'S 1 4 0 H*3V^T. iy=nuinn 

ord («f^<7)Ma7&^'l^t>ofc) ij^^ii^^m^-^ho :^t->;/7'S 1 4 0 C0^flJBff;&5Y 
ESC0t#^t;{i. -^ty^l^-^^^fett^o :^7^yrS l 4 0 <Dm^i)^^ O (Dt^'k 

5:^Oo 7.y^ yfS 1 4 2 tlioV^T. iy=numnord-l ({^^<?)*f^) T^ab^j^^SjO^ 
^mm-^ho 7sy- y':rs\A 2(r>m^'h^^0<Dm^\,Z{t^ 1 4 4 I'it 

tfo ;^x-7ys 1 4 2 (T^flJIIf/^'^^NOOi^^tCfi. 1 4 3 tlib^v^T. 

- K#-^idst) tl^i^o wluHC (7 - K#-^isrs:^-'^) tLT < ^ x*- • 

[0 0 8 1 1 

Xx-y^S 1 4 4-e(i, wlu2<7)T-'-^^j£o/cMif«0>i^vv3 >*ft»1-'l>o 
ibp=mod(ibp-l+numnord-l,numnord)+K ncptr=nbe+(ibp-l)*len+l ( 

[0 0 8 21 

t wlu2^5l§:iS-ro ^x <y ys 1 4 6 tCiSV^T. iy=iy+2^. 2^mx.X. 
;^-r y T'S 1 3 3 tCM'Bo 
[0 0 8 31 

112 4(i. •iryj\^~^ypim<Dyn-X'$)ho 

1 5 0lCi3i/^T. A (k. n/numnord) , wlul (lenblks. ibiksm 
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aero) 'i) L < ti^ wlu2 (lenblks, iblksmacro) ^wlux(lenblks, iblksmacro) 
{Z^if^o O^O^^mtx^i'^'bit^^fj 1^7Cg<7)^$&f]ASncptr=^^oTiE:^ff^(7) 

TiRimW^ir^o is2d=i*iblksunit+l^ ie2d=n/nuInnord^ len=ie2d-is2d+K isl 
d=ncptr. ield=nptr+len-l ( i {i-^f r^V-^ > p L Uco^ i9 ^ Llfc) . A (isld 
:ield, is2d:ie2d) =A (isld:ield, is2d: ie2d)-wlu(iblksmacro+l: iblksmacro+ 
len, 1: iblksmacro) X A (isld-iblksmacro: isld-1, is2d: ie2d) (^1) t^t^o 
[0 0 8 4] 

y^x y-7°s 1 5 1 tcisv^-c. 3k^\\ifzmm-^ ^ :^ y V^<D^1jU'^^i6XW '0 

_htf" numroot=int(sqrt(nuinthrd))> ^) L> sqrt (numthrd)-nuinroot:6^ 0 "C^fc 
U^^ib. numroot=numroot+li:1-'2>o «1«1T% i n Ht/h^j^MTW^t^X, s 
q r t ^yj^X-^^o :^7--yy°S 15 2 tC:feu>-C. ml=numroot. m2=nuraroot 
^ mx=ml i: i" <i) o X-r yT'S 1 5 3 twi3V>T> ml=mx> inx=inx-K imn=inxXm2 1 -i) 

0 T.T-y-i'S 1 5 4t3i3V>T. mm<numthrd-e^>;5;6^S7&^^r*IJ»f1"^o :^7-y-fS 

1 5 4<7)^flJWf75^'NO<7)^'g-lc{i, >?>-r -y/S 1 5 3 l^M^o X-r yT'S 1 5 4 <7) 
^iJifT&^'YE Scoi^-^tCti, 7.x-;'yS 1 5 5 tl^SV^T. Mffl" ^ fl^t ^ 1 >X7C B 

l^:g--r'i)o 2^7Cg^m2^53^L-C, m 1 X m 2 ^@<7)^§f^tc-t'2>o -^O^ 
numthrdli^^^Xl^-y KtCtiJ»5^-rr. (^1) O^SgP^^-^ilfe^iJl^ft^-r 
ho (1. l)^ (1, 2). ••• (l,m2). (2, 1) ••••<h2^ 

[0 0 8 5] 

>^x>y7'S156t343i/>T. ml*m2-numthrd>0;&^S>i^=^¥'JBIf-f'2>o Xx-v/S 
1 5 6<7)fiJBf7&^'YE S<7)Jt^>^tc{i, Xx ^T^S 1 5 8 Hjitfo J^t^ -^^-ZS 1 5 6 

<7)ft^^<7)^T> l^f<7)*f^75^ibml*m2-numthrd1@;6^'Mff$ti,-f t^^oTv^^o -::<7) 
^Eff^^lt^LT 1 o<7)^g?l^^#^, 2 ^xg 1^ y Kifcnumthrdt?55'#J LT ( 
^1) 0*tSg|5^^3fe?ijHlt^-r;2>o -eLT^ >?.xy7°S 1 5 8t;i3V^-c. /n''; 
TlBjA! (^l^yKW) ^i:o-C. -9" - ^ > felt's o 
[0 0 8 6] 

112 5{i. •9--7';v-^>fblu<7)-7n--e^;|>o 
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^.T-y-fS 1 6 OtCiSV^T. A (k. n/numnord) . wlul (iblksmacro. ibl 
ksmacro) ^ bufs (iblksmacro^ iblksunit) > bufd (iblksmacro, iblksunit) 

^^li^-e^ttr. - K(7)ipgiblksunit<7)ft|^coyn y ^ ^numnordtS^^fe/c ^ 
- K-cmil LT#o J: ^ HflJffl/f^J£^^^#y - Ktcii>i)o - Vi)"^ 
iblksmacroX iblksmacro<7)^n y ^ LT^o fjf^, ~ KT'lsI D^T^iJtC 

*tLTLU^Ii?^^f do LU^^jllT&^'^TL;^^. ^y- K ICSES ^tLTv^ytigp^ 

[0 0 8 7] 

^T'yfS 1 6 1 "Cfi^ nbase=n- ibl ksmacro ^ ibs=nbase+K ibe=n^ len=iblk 
smacro^ nbase2d=(i-l)*iblksunits ibs2d=n/numnord-iblksunit+l> ibe2d=n/nu 
mnordi: L^ iiigT^'— :5!'iSL{ilensend=iblksmacro*iblksunit ^ iy=l<i:"t'-?)o 
[0 0 8 8] 

;^-r yT'S 1 6 2 tcjsv^Tti. 7 y^cT) n fcf - ^^f d o i-^t>*>. bufdd 
:iblksmacro, 1: iblksunit)— A (ibsribe, ibs2d:ibe2d) h-t"^o ^-r yT'S 1 6 
StCioV^Tti. iy>numnord75^S>0^*flB/fi-^o y^x y^'S 1 6 3 OfiJ8lf7&^~Y E S 

t>%^ idst=mod(nonord-l+iy-l, numnord)+K isrs=mod(nonord-l+nuinnord-iy+l, 

numnord)+l<hi->&o T.y- yy'S 1 6 ^7- KTii^fti-^o bufdC0I*I# 

^idst#g(50y - Ktcii|.i>o X-f y^^S 1 6 6t;i3V^T{i. bufsilr"*- ^ 

y-fS 1 6 8tli3\/^-r. wlulc7)HCDfi«Cisrs#g<oy- Vi)^h^fzr- 9 
^i"^o icp2ds=(isrs-l)*iblksunit+K icp2de=icp2ds+iblksunit-l> wlu(l:ib 
Iksmacro, icp2ds: icp2de)— bufs(l : iblksunit, 1: iblksunit) h'fho ^T-yy'S 
16 9i;^v^T. iy=iy+l^L. -r y/S 1 6 3 tC^'5>o 
[0 0 8 9] 

X7"y7°s 1 7 o-e{±. rmm^kt 19. ;=^7" 7 7°s 1 7 1 -eti. wiuioji 

-eiblksmacroX iblksmacro<7) L U^M^^y ~ K-CS^ LT^t9 o tf ^^<7)1t^ 
ti^ ipt3;^^-t-^o LU^jST&^ll-TLT^c'b. g ^ - K^^ftm<7)-rn y n k: 
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— /"i y ^ir^o 't"'&t5*>> is=(nonord-l)*iblksunit+l, ie=is+iblksunit-l^ A 
(ibsiibe, ibs2d:ibe2d)^wlul(l:iblksniacro, is:ie) ^ LT> 

(#151) m^(7)-fu^y^t^^') ^^ism^<DJ- Vi^^ y hy-i^t? 

^tiny^u y ^ - Kt;%fttc@etti-.i> %%nfu y^^i >k7QB 
K-e^^^iJlc L U55-|lli-'i> L H^m^ T'-y-i't, 
[0 0 9 0] 

[0 0 9 1 1 

(#f&4) BufEife^iJft^^tix SMP (SymmetricMultiProcessor) 
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5 ) y'u-ty^t ^^') y - K ^ y h 7 - ^ -e 

L fzmnmm. it ^ ^oM^s-e^ o r . 
yxmi] KLu ^m-t ^L\j ^m^m. t . 

[0 0 9 3] 

y - K^tHt-r -7 t3#J «9#tt 'btL/itf^iJ<7)g|5:^<7)^ij7'n y^(Dim^^ 

m\m^5^^ttt>tzii<D^iMmz\^x^mthfzib\z^ #y-Ktc— o-r 

X6^m-r h^t\z^^xmhKhy^u y ^ "kmW^fD y - ki;. ^tcii&^iJilffi L 
-C-o-roKSi-.|,^2 OKgx-r -y y<!:. 
ig*t:^-7*n y tgea^tLfc-rn y i^^, #y - KrelT^a^ L^T^^'ib. #y- 
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1 # § 55- ^ IS-a- L 7~ n y ^ tcM L T :^ fl^^ h f u y ^ (n> Tf « 

mnny^ y ^ - Kn%ftnBEg-r'2> t^^^ mmr^ y^n >X7tE 

x-'^-r-Dum.-f^m2mim.^7-y y°t. 
mnnf'ci yi7twM.^i\^fzy'U 'y'^^, - vrnxmrnL^^-h^h. 
Yxm^i L u ^^m-f ^m:^ f-yft. 

-ft. 
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[0 0 9 7] 

t X\ MM y^ii -^i6T S M P (Dskmit^^X # T S M P X(7)^m^mx<7)it 

[mi] 

[1112] 
[0 3] 
[[114] 

m (-^(D 1 ) x$>^o 
[ins] 

m (-e<7)2) -c^^o 
me] 

imi] 

nmm^j:Lu^m<D^mim.m-r^mx$>^o 

[0 8] . 
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[119] 
[HI 0] 
[Hill 
[H 1 2] 

mi 3] 
mi 4] 

[HIS] 

[Hie] 

mi7] 

^§tm<^mM^^w.<7)yu~^-^~Y (^<D6) x^^o 
mi 8] 

:$^§^m(ommMm<7)yu-^^~h (^(d?) x^^o 

[HI 9] 

m2 0] 

[(112 1] 

^^m<^nmfm<^y^-'^^- h (^(dio) x$)ho 

[1112 2] 

^^m(DmMfm<oyu-^^~}> {^<di d x^^o 
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[1112 3] 
[12124] 

[E12 6] 

1 0 tiSlS-g-ffl (z^-;^) 

1 1-1-1 l-n 

12- l~12-m ^^yvj. 

13- 1-13 -m 

14 x-^iijfffl/N- K (DTU) 
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K3 
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7, 



B2 
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24 
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B6 



B3 



B7 



II 



y-Ki 
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/-K2 
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imi 0] 




/^y3'T— B 
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imi 1] 



J... ; .. .,-1 L 
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imi 2] 



numnord <!:L'Crv=iblksunit x numnord X m (mli^> — K'^0-^(i!^P*>^ 
(k.n/numnorcl) (k>=n) feJ:t;^t©X+lBSx.©HS*ft<|l|1-Sip(n) $?IS8t 



S20 



I 



nonordl=^D-fe::<.«#( 1 K»)$SS. 



S2t 



\ 



S22 



tblksmacro=ibIksunit x numnord yu^y<7^mO^^ 

loof>=n/(iblksunit x numthrd)-1 (SiygLlHiSS) 

i=l . lenbufmax=(n-iblksmacro)/nunnnord+ibB^smacro 



T 



wluKlenbufmaxjblksmacro) wlu2(lenbufrnax,iblksmacro) 
bufsdenbufmaxjblksunit) bufddrnbufhnaxJblksunit) 

lenbuf£tt-»Crjes®/d:*#*«lt«5» 



S24 



S25 




S26 



T 



lenblk$=(n-i x jblksfnacro)/nuinnord+iblksnriacro 



Kf::fe5>*Biblksijnit <0\ 
ibIksmacroCO:/ □ vO^Se^fk ^ 



S28 



T 



S29- 



^^interhJ ^I^Utli 

is=(i-1 )*iblksmacro+ 1 . 
ie^HHblksmacro tL"Cip(isic)Iw 
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S30^ 
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S37 
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S35> 



S36 
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fblu^l^^^?g 
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*^BJ(7>||M^»(D:7P-5^-V-h(-tCD3) 



S45 



A(k.n/numnord), wlu 1 (lenblksjblksmacro) 
bufsdenblks.ibiksunit) bufd(lenblks,iblksunit) 
^lliST'SttT. S-y-KOi ||a0i|giblksunit (DZf Uv^^ 

numnord fil^^WLfct^i^^ft^^O'V^^an^fctt^^y—Kl: 



S46 



nbase=(i-l )>»=iblksmacro, 

ibs=nbase^1 ,ibe=nbase+iblksmacro 
len=(n-ibe)/numnord ,nbase2ct=(i-l )*iblksunit 
ibs2d=nbase2d+1.ibe2d=ibs2d+iblksumt 
SHItt — ^|K(^lensend=(len-i-iblksmacro)*lblksunit 




S49 



idst=mod{nonord-l +iy- 1 ,numnord)+1 
•srs=mod(nonord-1 +numnord-iy+1 ,numnord)+1 



S50 



ft/— K-Ce»lwf8yoLxTL^«ffiiblksunit C0»^^:;^O-:/^7 
SP^^t. ^■©TgB5J^CD::?0'V^(7)i::)^7£Bi&numnord T-^^ 

bufdCl :iblksmacro.1 :ibiksunit)'^A(ibs:ibeJbs2d:ibe2d) 
icps=ib€+(idst-1 )*len^^l .icpe=isps+len-l 
bufdCiblksmacro+l :len+iblksmacro,1 :iblksunit) *— 
A(icps:icpejbs2d:ibe2d) 



S5I 



\ 



SISmt«.(^>'-KT-2Hg>bufd(3[>rtS$idstSBO!>y- 

Kicai-5bufsizsfi-r€>o 



0 
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[lai 5] 
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S52 




S53- 



S54 

VI 



wiui <on!$i^m.iz^^-t^o 

isrs# S CDy--K3&^C>SJ+fcT-4'^*»I 

icp2ds=(isrs-1 )*iblksunit+1 . 

icp2de=icp2ds+iblksunit-1 

wlu1 (1 :ien+iblksmacro,icp2ds:1 cp2de}<— 

bufed :len+iblksunit1 :iblksunrt) 
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S55 



iy=iy+l 



C 



return 
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\Y A(k,n/numnord), wlul(lenblks.iblksmacro) \ 
\, wlumicro(ncash) (wlumioro$'L2+-\"yi/^a.(Z)::^#$ J 




S62 



\ 



iblksnnicro=nwiclthmicro 

P'^'^A^I&Srt^F^lTt^^wiudenmacro.iblksmacro)© 
wlu(istmicro:lenmacro,istmicro:iblksmicro+iblksmicro-1 ) 0) 

wlu(istmicro:istmicro+iblksmicro--1 , istmicro: 
istmicro+iblksmicro-1 )^d\agil^^o 
irest=i5tmicro+iblksmicro 

wlu(irest:lenmacrojstmlcro:istmicro+iblksmicro-1 )^1 <)i 

U^>Kc:i:(^>l|iawlumicrolZZie— -^So 
lenmicro=(lenmacro-~irest+numthrd)/numthrd 
wlumicrodenmicro+iblksmicrojblksmicro ) IwUtz — o 
lenblksmicro-lenmicro+iblksmicro 



S63- 



cro^WtB-To 
( wlumicrodenmicro 
+iblksmicrojblksmic 

ro )^sitai-ro ) 
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wlumicrofC^iljLrt^fcSP5J$. diag ^$3^1*100 
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nwldthmicro2=nwidthmicro/2 
istmi cro2=i stm icro'i'nwidthm i cro2 
nwidthnriicro3=nwidthmicro-nwidthmicro2 



S67 



nwidthmicro2=nwidthmtcro/3 
istrnicro2=istmicro+nwidthmicro2 

nwidthmicro3=nwidthmlcro-nwidthmicro2 



S68- 



istimicro I4'?:<7>$ 

nwidthmicro tLT 
nwidthmicro2£^SiLT 



S69 



wl u(istm icroti stmacro-t-n widthmi cro2- 1 . 
istmacro+nwidthmicro2:istmacro+nwidthmicro- 1 )ff> 

5?) u.itlCwlu(tstmicro:istmacro+nwidthmicro2-1 , 
1 stm i c ro:i stmacro+nwidth m i c ro2- "F^ ^1 tT?*J © 



S70 



\ 



T 



wlu(istmicro2:lenmacro. 
istmicro2:istnnicro2-hiwidthmicro3-*t 
wlu{istmicro2:lenmacro.istmlcro!istmicro2-1 ) x 
wlu(istmicro:istmacro+nwidthmicro2-l, 
istmacro+nwidthmicro2:istmacro+nwidthmicro*l)^ 



istimicro iLT 
istmicro2 ^ nwidthmicro 
irLTriwidthmicroSS-iS 



return 
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[mis] 



A(k.?i/numnord) wlutOendks.itlKsmacro) 
wlumicro(leniblksmicro.tblksiTiicro)(wlurnicro^ 

:|:jU-^>t:wlumicrol=»»$ntc«»a)LU»«f$?T/dr5c 



/ 



S75 




S79- 



S80 



T 




S85 



S83 



\ 



T 



S66 
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[019] 

0 




S89 



nwidth2=nwidth/2 
ist2=ist+nwicH:h2 
n wi dth3=n width-n width2 



nwi dth 2=n width / 3 
ist2=ist+nwidth2 
nwidth3=nwidth'*nwidth2 



S90 



S91- 



nwidthtL 
"Cnwidth2 ^ig 

^>LUmicro 



S92 



\ 



wluTnicro(istmicro:istmacro+nwidth2-l . 
Istmi cro+n wi dth2:i stmicro't-nwi dthmicro-l )© 

wlumicro(istmicro:istmacro+nwidth2- 

1 ,istmicro:istmacr(rHiwidth2-1 )(f>T=.fkff9i\ 



S93 



I 



wlumicro(istniicro2:lenmacro,istfnicro2:istmicr 
o2+nwidthmicro3-1 
wtumicro(istmicro2:lenmacro. 
istmicro:istmicro2-1 ) x wlumicro 
(istmicro:istmacro+nwidth2-l , ist+nwidth2: 
ist+nwidthmicro-1) $§|l^"CMSfft*S 
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[1112 0] 



SI 00 



■ ■ 

1^ / A(k,n/numnordj. vuLiKlenbiks.iblksmacro) bufsdenblks.ibtksunit) 

\ / bufdClenblks.iblksumO^^li&trSftT. #y--K©i Sia>«|iSiblksun!t (Oyay<7^ 
^ numnord <l^tefc^»f^O)5^^33^a»>^f75»l§e^iblksmacro Xiblksmacro^fcU 




stoi 



\ 



ibs=nbase'*-1 jbe=nbase'i-iblksmacro 
len=(n-ibe).''numnord ,nbase2d=(iH >iblksunit 
ibs2d=nbase2d^1 .ibe2d=ibs2d«iblksunit 
^ tex — ^ S!( lensend=( len+iblksmacro)*iblksunit 




idst=rno<l(nonord-t +»y-1.numnord)+1 
isrs=mod(oonord- 1 +numnord-iy+ 1 .numnord)+ 1 



S105 



\ 



icp2ds=(id5t-1 )*ibM(5unit+1 . jcp2de=icp2ds+iblk5unit-1 
bufd(1:lerr»'ibtksunit.1:iblksunit) *- wlu1(1:lerH-ibtksmacro.icp2dsJcp2de) 
1 3fe5cB $;^U-vK»T-»WUT&:7.U-^Kr*3fe5»j3lf--r 5. 



SI 



bufda)|^$$idstSS 43r>/~Kl=]£^ 



St07 



SlOi 



I 



S109 



A(ibs:ibe.ibs2d:ibe2d>~buf5(1 riblksmacro.l :iblksuni0 icpssibe-Kisrs-l )*len+1 . 
icpe=i5ps+len-l 

A(icps:icpe.fbs2d:ibe2d> — bufsCiblksmacro^l rlen-t-iblksmacro.l :iblksunit) 



S11 



c 



return 
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2 1] 



SI 15 / A(Kn/numno*d). wluKlenblks.iblksmaco) 

wlu 1(1 :iblksmacro, 1 liblksmacro) iZitLV^M^tltz 

nbdiag"(f-0*iblksmacro (WtmiS^HLTtO-^^Zf 
ip(nbdiag+1 :nbdiae*-iblksmacro ) ICftlft^tuTL 



S11 



T 



nbase=i*ib!ksu nil 

irows=nbase+1 .lrowe=n/numnord 
ten-(irowe~irows+1 )/nurrithrd . 
is=nbase+(nothrd- 1 )»lef»*l , 
ie=fnin(irowe js+len-1 > 



8117 



sua 




nbdiag=(i- 1 )*iblksmacroj=nbdiag- 



Ati.ix) ^A(tp(j)Jx) ^Xtim^L^ 



SI23 



S119 




ix=ix+1 



S124 



3125-^ 
8126^ 



A(nbdiag+1 :nbdiag+ibfksmacro.is:ie) ^ 
TRUwiuld :iblksniacro, I liblksmacro))- 1 x 
A(nbdiag+l:nbdiae*tbIksmacro.is:i8) $ 



8127^ 



c 



return 
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wluKlenblks.iblksmacro) wlu2(ienblks,iblksinacro} 



S131 



SI 32 



nbase={i-1 )*ibIksmacro. 

ibs=nbase+1 .ibe=nbase-4-iblk&mdcro 
len=(n-ibe)/numnord ,nbase2d=(i-1 )*iblksunit 
ib52d=nbase2d+l4be2d=ibs2d+ibiksunit 
n2d=n/numnord 
lensend=len+iblksmacro 
S^fg-X — ^ IS 14 nwlen=lensend*iblksmacro 



iy=l(«]»!m^iSS) 
idst=mod(nonord.numncrd>+l< iSiyife/— KS-^(B?/— K) ) 
isrs=mod(nonord-1 +numnord-l ^numnordHl ( 7C-/-^KS%) 

ibp=idst 



SI 33 



5135- 




S126 



SI 37 



8138 



\ 



ibp=nnod(ibp-l +numnord-1 ♦numnord)+1 
ncptr=nbe4<ibp-1)*len+l ( 1 ;ji^5cBCDSaj!fe{ifi) 
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S142 



SI43 




•f-isrs*^c,)5ie,tt-c<*7^-^i$*ft«-rs, sL&m 

■T— ^felinwlen „ 



S144 



ibp=mod(ibpH +nurnnord— 1 ^umnord)+l 
ncptr=nbe-Kibp-1)*len+1 ( 1 JfejcBtDBBtefiB) 






iy=iy+2 (2 SjtB!3.S) 







c 



return 
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A(k,n/nuninord), wlul (lenbIks.iblksnriacro)tiL< H 
w]u2(lenblks.ibiksmacro)$wlux(lenblks,iblksmacro ) iZ^i'f^o 
l>¥l/tt5Lti)i^?>a**tfc1*7cBCZ>BB*Safincptr ^ttoXiE 

is2d=i*iblksunit^-1 .le2d=n/numnord Jen=ie2dHs2d+l 
is 1 d= n c ptr. i e 1 d=n ptr+i e n- 1 
(i li-9-:7Jl^-5^>pLU®ISyiSLa) 
A(ist d:iel d.is2d:ie2d>=A(is1 driel d,is2d:ie2dh 
wlux(iblksmacro+l:lblksmacro+len,1 :iblksmacro} x 
A(i$ld'-iblksmacro:isld-1.is2d:ie2d) (jECI ) 




S151 
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I 



numroot = int(sqrt(numthrd)) 
iKsqrt(numthrdV-numroot.ne.O) numroot=numroot +1 



S152- 
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SI 54 
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m 1 =numrootm2=numroot 

nix=m1 














ml=mx 
mx=mx-1 
mm=mx Xm2 










-<C^7 mm<nufnthrd 



S-;?.U'>Rcsiiy^rr(Si)oJtji£;gi5^^^at5>jJctt-e-r 

So (1 .1 >.(! .2)....(1 .in2).(Z1 )...t2J»[7t S ®* iSllZflSSJCX 
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S157 
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ml*m2-numthrd ^t^fSit^'flt^-oXl'^ho 

«numthrd r^fflUTCsci )<D«i5»»^afe5ijf=stir-r 
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Alk.a'^nunrinord), wlut (iblksmacr u.iblksma^vD) 
bufs(iblksmacro.ib!ksunft) bufdOblksmacro^iblksunit) 

^y—K i<iblksmacro x ib!ksmacroCD:?P*:^i?$'SltUrlto/iS. 
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nbase^n'tblksmacro. 
ibs=nbase-»-l .ibe=n 
len=ibtksni<tcro .ntiase2d=(i~1>^btksunH. 
ib&2d=n/nuninord-iblk&unit+l,ibe2d=n/numnord 
vk^T^ — $SSt(ilensend=iblk&macro*ihi<sunrt 

iy=i 
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id8t=mod(nonord~'1 -Hy~1 .numnordH'l 
isrs=nnodCnonord~-1 +numn(>rdSy+1 ,numnord)+1 
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icp2ds=(isrs- 1 )*iblKsunit+l . 
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